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The ultrafast relaxation and recombination dynamics of photogenerated electrons and holes in 
epitaxial graphene are studied using optical-pump Terahertz-probe spectroscopy. The conductivity 
in graphene at Terahertz frequencies depends on the carrier concentration as well as the carrier 
distribution in energy. Time-resolved studies of the conductivity can therefore be used to probe the 
dynamics associated with carrier intraband relaxation and interband recombination. We report the 
electron-hole recombination times in epitaxial graphene for the first time. Our results show that 
carrier cooling occurs on sub-picosecond time scales and that interband recombination times are 
carrier density dependent. 
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Graphene is a 2D lattice of carbon atoms arranged 
in a honeycomb crystal structure with a zero (or near- 
zero) bandgap and a linear energy-momentum disper- 
sion relation for both electrons and holes P, 0]- The 
unique electronic and optical properties of graphene 
make it a promising material for the development of 
high-speed electron devices, including field-effect tran- 
sistors, pn-diodes. Terahertz oscillators, and electronic 
and optical sensors [3, H, 0, S, @, 0] ■ The realization of 
graphene-based devices requires understanding the non- 
equilibrium carrier dynamics as well as the rate at which 
electron-hole recombination occurs. 

Measurements of the ultrafast intraband relaxation dy- 
namics of photogenerated electrons and holes in epitaxial 
graphene using both degenerate Q| and non-degenerate 
[9| optical-pump optical-probe spectroscopy have been 
previously reported. Similar measurements for exfoliated 
graphene mono- and multi-layers have also been carried 
out [lo| . These measurements were sensitive to the in- 
terband conductivity of graphene and probed the time 
evolution of the carrier occupation at specific energies 
in the bands. Consequently, they were not able to di- 
rectly measure the time scales associated with carrier re- 
combination. At room temperature, the optical response 
of graphene in the THz frequency range is described by 
the intraband conductivity - the free carrier response - 
which depends not only on the total carrier concentra- 
tion but also on the carrier distribution in energy 11 1. 
Therefore, THz radiation can be used to study the car- 
rier relaxation and recombination dynamics in graphene. 
In this paper, we present results obtained from optical- 
pump THz-probe spectroscopy of epitaxial graphene in 
which the time-dependent conductivity of graphene that 
has been excited with an optical pump pulse is probed 
with a few-cycle THz pulse. We observe cooling of the 
photogenerated carrier distribution as well as electron- 
hole recombination in graphene in real time. Our results 
indicate that the recombination times in graphene de- 
pend on the carrier density and material disorder. 

The epitaxial graphene samples used in this work were 



grown on the carbon face of semi-insulating 6-ff-SiC 
wafers using techniques that have been reported previ- 
ously [13]. As discussed in @, [ll|, X-ray photoemission, 
Raman, and optical/IR/THz transmission spectroscopy 
were used to characterize each sample to determine the 
number of carbon atom layers and the carrier momen- 
tum relaxation time. In this work, sample B consists 
of ^^12 carbon atom layers with a momentum relaxation 
time r ~ 20 fs and sample C consists of ~19 layers with 
T ~ 4 fs Uj. Raman spectroscopy (at 488 nm) of the 
samples revealed a G peak close to 1580 cm~^ and a D 
peak close to 1350 cm~^. The D peak is forbidden in per- 
fect graphene layers since it requires an elastic scattering 
process, which is made possible by disorder, to satisfy 
momentum conservation l^, 3] ■ The presence of the D 
peak therefore indicates the presence of disorder in the 
samples. The average ratio of the intensity of the Raman 
G and D peaks, Iq/Iy^, was 17 for sample B and 2 for 
sample C. This indicates significantly more disorder in 
the latter sample [l^ . Results from optical-pump 
optical-probe spectroscopy of sample B were reported 
earlier [8|. 

Ultrafast 90 fs pulses from a 81 MHz Ti:Sapphire laser 
with a center wavelength of 780 nm (1.6 eV) were used 
to simultaneously pump the THz Time-Domain Spec- 
trometer (THz-TDS) and the epitaxial graphene sam- 
ples. The THz-TDS used a photoconductive emitter and 



an electro-optic detector [15|, |16[, and generated few- 



cycle THz pulses with a 0.3-3.0 THz bandwidth and a 
power SNR near 10^. The epitaxial graphene samples 
were placed at the focus of a 90° off-axis parabolic mir- 
ror and excited with 1 nJ - 16 nJ optical pulses focused 
to FWHM spot sizes of 350 /zm - 400 /xm. The time- 
dependent transmission of the THz probe pulse through 
the optically excited sample was measured by delaying 
the optical pump pulse with respect to the THz probe. 
The pump and probe signals were mechanically chopped 
at 256 Hz and 179.2 Hz, respectively, and the change in 
the transmission was measured using a lock-in amplifier 
referenced to the sum of these frequencies. 
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The complex amplitude transmission of the THz pulse 
through the graphene layers and the SiC substrate (nor- 
malized to the transmission through the SiC substrate) 
is, 



t = 



1 



1 + Nr]o (Tintra/(1 + "Sic) 



(1) 



N is the number of carbon atom layers, r/o is the free- 
space impedance, and nsic is the refractive index of SiC 
and equals ~2.55. The transmission is related to the com- 
plex intraband conductivity of graphene, Uintra, which is 
given by 



f^intra 



i^-— / [/ (6 - ep„ T) + f{e + epv, T)] de 

^ + i/T Jo 

(2) 

In this expression, t is the momentum relaxation time, 
/(e — eF,r) is the Fermi-Dirac distribution with carrier 
temperature T, and epc and epv are the Fermi energies 
of conduction band and valence band electrons, respec- 
tively. The electron and hole densities in graphene are. 
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Equations [2] and [3] show that, unlike in semiconductors 
with parabolic energy bands, the intraband conductivity 
of graphene depends on the carrier distribution in addi- 
tion to the carrier density. In graphene, a hotter carrier 
distribution corresponds to a smaller intraband conduc- 
tivity than a cooler distribution with the same carrier 
density. In intrinsic photoexcited graphene layers, n — p 
and epc = —epv In epitaxial graphene, the first one or 
two carbon atom layers adjacent to the SiC substrate are 
expected to have a large carrier density [^, How- 
ever, the change in the THz intraband conductivity, as 
measured in our experiments, is dominated by the large 
number of intrinsic layers. The contribution of the doped 
layers to our measurements has therefore been ignored in 
the discussion that follows. 

Figure [T] shows the measured THz pulse transmitted 
through graphene sample B without and with (scaled for 
clarity) optical excitation for a probe delay of 1 ps. In 
this work, we measured the change in the transmission of 
the peak of the THz probe pulse caused by optical exci- 
tation. Because of the small momentum relaxation time 
of samples B and C, the real part of the intraband con- 
ductivity is much larger than the imaginary part and is 
nearly dispersionless in the low THz frequency range. In 
addition, the first order contribution from the change in 
the imaginary part of the conductivity to the change in 
the pulse transmission is exactly zero at the peak of the 
THz pulse. Consequently, near the peak, optical exci- 
tation primarily affects the amplitude of the THz probe 
pulse and does not distort the pulse shape or shift the 
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FIG. 1: Measured THz pulses transmitted through the epitax- 
ial graphene sample B without (grey) and with (black, scaled) 
an optical pump pulse preceding the peak of the THz pulse by 
1 ps. The the intraband conductivity of both samples B and 
C is almost entirely real and is nearly dispersionless at low 
THz frequencies. Consequently, the amplitude of the peak of 
the THz probe pulse is altered in the presence of the opti- 
cal excitation but the pulse shape is not distorted nor is it 
shifted in time. Ringing after the pulse is due to absorption 
from water vapor. 
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FIG. 2: (a) The measured change in the real part of the 
complex amplitude transmission (grey) and the theoretical 
fit (black) for pump pulse energy of 14.8 nj (sample B). The 
transmission decreases rapidly until 0.75 ps - 1 ps and then 
increases slowly from 1 ps - 15 ps. The inlay shows the initial 
rise in the transmission, which occurs on a timescale longer 
than the experimental resolution of 130 fs. (b) The same 
data plotted on a logarithmic scale. The slow increase in the 
transmission does not follow an exponential curve. 
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FIG. 3: A schematic of the hkely processes by which optically- 
excited, non-equilibrium electron and hole distributions ap- 
proach equilibrium. After excitation, the distribution rapidly 
thermalizes and cools. The hot thermally distributed carriers 
are then cooled further due to intraband phonon scattering. 
Finally, electrons and holes recombine until the equilibrium 
distribution is restored. 



pulse in time. Figure [Ifa) displays the experimentally 
measured time-dependent change in the real part of the 
normalized transmission of sample B resulting from ex- 
citation by a 14.8 nJ optical pulse. The time evolution 
of the transmission exhibits two distinct features: (1) an 
initial rapid decrease until 0.75 ps - 1 ps (left of dashed 
vertical line in inlay) followed by (2) a slow increase from 
1 ps - 15 ps. The initial rapid decrease, shown in the in- 
lay in Fig. [21 is within the (estimated) 130 fs temporal 
resolution of our measurements and is consistent with 
an increase in conductivity due to cooling of the ther- 
malized carrier distribution after photoexcitation. The 
slow increase in the transmission is due to the decrease 
in conductivity resulting from electron-hole recombina- 
tion and does not follow a simple exponential dependence 
(Fig. Etb)). This is indicative of density-dependent car- 
rier recombination times. 

The likely processes by which non-equilibrium photo- 
generated electron and hole distributions return to equi- 
librium are shown pictorially in Fig. [3] Immediately fol- 
lowing optical excitation, electrons (holes) are in a non- 
thermal distribution with a large population near 0.8 (- 
0.8) eV from the band edge. Within 10 fs - 150 fs, this 
distribution thermalizes and cools to form a hot Fermi- 
Dirac distribution Intraband phonon scattering fur- 
ther cools the thermalized carriers between 150 fs and 
1 ps. Beyond ~1 ps, electron-hole recombination is the 
dominant processes affecting the intraband conductivity 
and the THz transmission. Intraband relaxation times 
in the 0.4 ps - 1.7 ps range were measured by the au- 
thors previously using optical-pump and optical-probe 
spectroscopy of epitaxial graphene Q. The intraband 
relaxation times measured in this work are slightly dif- 
ferent. This small difference can be attributed to the 
fact that in Q the time-dependent optical transmission 
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FIG. 4: The theoretically predicted and experimentally mea- 
sured changes in the real part of the transmission for optical- 
pump pulse energies of (a) 9.9 nJ, (b) 7.4 nJ, and (c) 4.9 nJ 
(sample B). The model captures both the fast decrease and 
the slow increase of the transmission for different pump pulse 
energies. 



depended on the carrier occupation at a specific energy 
in the high energy tail of the distribution, whereas in 
the present case the time-dependent THz transmission is 
sensitive to the changes in the distribution at all energies. 

Simple coupled rate equations can be used to model the 
carrier dynamics after optical excitation and subsequent 
thermalization. These equations model the cooling of 
the carrier temperature T towards the equilibrium tem- 
perature with a phenomenological time constant and 
the decay of the carrier population due to recombination. 
The mechanisms responsible for electron- hole recombina- 
tion in graphene could include plasmon emission, phonon 
emission, and Auger scattering. In general, recombina- 
tion rates due to these processes are expected to have 
non-trivial carrier density dependencies [3, [20l . 



21, m 



Additionally, carrier generation rates cannot be ignored 
in the analysis 2l|. We assume that the net recombina- 
tion rate R{n) can be expressed as a quadratic function 
) , where it is assumed that n = p 



-eq 



of n, R{n) = B{n^ 
and rioq is the carrier density at thermal equilibrium. B 
is a fitting parameter. From the time-dependent carrier 
temperature T and density n, the conductivity Cintra and 
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FIG. 5: The measured and the predicted peak change in the 
real part of the complex THz transmission for different pulse 
energies (sample B). 



the THz transmission can be determined. 

The agreement between the model (black) and the 
measured data (grey) for the THz probe transmission 
transient is also shown in Fig. [21 For a pump pulse en- 
ergy of 14.8 nJ, the photoexcited carrier density was es- 
timated to be ~ 5 X 10^^ cm~^. To best fit the measured 
data, the carrier temperature immediately after thermal- 
ization was assumed to be 1250 K, and the values of Tc 
and B were found to be 200 fs and 1.8 cm^ s~^, respec- 
tively. To explore the carrier density dependence of the 
measured transients, experiments were conducted with 
optical pump pulse energies varying from 1 nJ to 15 nJ. 
The model also exhibits an excellent fit to the data for dif- 
ferent pump pulse energies, as shown in Fig. [4] for pump 
pulse energies of 9.9, 7.4, and 4.9 nJ. In the model, the 
initial photoexcited carrier density and the initial carrier 
temperature were scaled linearly with the pump pulse en- 
ergy Q . It is important to emphasize here that the model 
captures the initial decrease in the transmission, the peak 
change in the transmission, and the slow increase in the 
transmission for different pulse energies without chang- 
ing the values of B ot Tc- Figure [5] shows good agreement 
between the measured and the calculated peak change 
(absolute minimum, see inlay) in the THz transmission 
for pump pulse energies varying from 1 nJ to 15 nJ. The 
extracted value of Tc agrees well with the calculated and 
measured inelastic carrier-phonon intraband scattering 
times in graphene [13] and in carbon nanotubes 18, 1^ . 
Using the extracted value of B, the calculated recombi- 



nation rates for carrier densities in the 10^^ — 10^^ cm~^ 
range are approximately 1.5-2.0 times larger than the 
theoretically predicted recombination rates due to Auger 
scattering in graphene 2l|. This difference can be at- 



expression. For example, in ^2)\ it has been shown that 
the Auger recombination rate in graphene can have a cu- 
bic dependence on the carrier density. The coefficient B 
can also be carrier density dependent 2l|. However, for 
the values of the carrier density relevant to our experi- 
ments the calculated auger recombination rates have an 
almost quadratic dependence on the carrier density 21 1 . 
In addition, within the resolution of our experiments, a 
constant value of B adequately fits the data. 

In and 11 1 it was shown that the intraband energy 
and momentum relaxation rates in epitaxial graphene de- 
pend on material disorder. Disorder could also affect the 
interband recombination rates since disorder can modify 
the density of states [23| and/or provide the additional 
momentum necessary to satisfy the conservation rules in 
recombination processes (23 . [25j . Optical-pump THz- 
probe spectroscopy was performed on epitaxial graphene 
sample C which, as discussed earlier, is significantly more 
disordered than sample B. The value of B for sample 
C was obtained using the methods discussed above and 
equals 3.1 cm^ s^^. The carrier recombination rate is 
approximately twice as fast in sample C than in sample 
B for the same carrier density. 

In conclusion, we studied the ultrafast relaxation and 
recombination dynamics of photogenerated carriers in 
epitaxial graphene using optical-pump THz-probe spec- 
troscopy. We measured the interband electron-hole re- 
combination rates in graphene for the first time. Our 
measurements show that the recombination times in 
graphene depend on the carrier density. 
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tributed to additional recombination mechanisms, such 
as plasmon and phonon emission. It should be mentioned 
that the net recombination rate R(n) can have a carrier 
density dependence different from the assumed quadratic 
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